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IN THE CLAIMS 

This is a complete and current listing of the claims, marked with status identifiers in 
parentheses. The following listing of claims will replace all prior versions and listings of claims in 
the application. 

1. (Currently Amended) A Method method for sampling a state space by iteratively generating 
states x*,t and their weighting factors p,t , wherein the index i is the iteration parameter and 

the index t distinguishes different states Xi,t generated by an iteration i, the method is 
consisting o f comprising: 

a first step for selecting an initial sampling distribution function p(x) , 
a fifth step for performing an analysis 
and an iteration procedure including 

a second step for generating Nj states xj,t by a numerical sampling algorithm and 

a fourth step for testing at least one criterion to decide whether to continue the iteration 

procedure or to stop the iteration procedure and to go to a fifth step in order to perform the 

analysis using the simulated data, wherein 

characterized In that the iteration procedure further includes 

a third step determining weighting factors pi,t for states Xi,t generated so far by using 
sampling distribution functions p. (x) determined so far and 

a fitting step for determining a sampling distribution function p • (x) for the next iteration by 
fitting Pj{x) to pi,tO(xi,t) for states xtj generated so far, wherein 0(xij) is a function, 
respectively a property, of the states xm . 

2. (Currently Amended) Method as claimed in claim 1, wherein the sampling distribution 
function pj(x) of at least one iteration is fitted such that it maximises maximizes an objective 
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function preferably defined as a function of local comparisons between the sampling 
distribution function and the product pi,tO(xij) . 

3. (Currently Amended) Method as claimed in claim l-ep-2, wherein the sampling distribution 
function pj(x) of at least one iteration is fitted such that the sampling distribution function 
pj(x) is large for at least one state Xij with a large product pjO(xi,t) , and tends to be small 
for states with a small product p,tO(xij) . 

4. (Currently Amended) Method as claimed in one of claims claim l-4e-3, wherein the 
sampling distribution function pj(x) of at least one iteration is a function with at least one 
constraint 

, preferably having an extreme value in tho region of at least ono selected state Xj, C onstr and having 
vanishing values at a distance from the at least one selected state x^constr , wherein the 
constraint is preferably a harmonic constraint with at least one constant -kconst F defining the 
strength of the constraint at the selected state ~xj ,constr 7 

5. (Currently Amended) Method as claimed in claim 4, wherein the sampling distribution 
function pj(x) of at least one iteration is the distribution function of the system with 
constraints or a multicanonical distribution function with constraints. 

6. (Currently Amended) Method as claimed in claim 1 or 2 , wherein the numerical sampling 
algorithm of at least one iteration generates correlated states xjj ; wh e rein the sampling 
distribution function pj(x) preferably has a maximum that biases sampling into a sub region 
of the state space, wherein tho sampling distribution function pj(x) is preferably fitted such 
that tho maximum is in a region where the product p t ,iO(xt,t ) has a maximum, and wherein 
the sampling preferably starts from a state close to the maximum of the sampling distribution 
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function - pj{x^. 

7. (Currently Amended) Method a& claimed in on e of claims 4 to 6 , wherein the fitting of 
pj(x) is done by selecting states x it for which the product p,tO{xi,t) has extreme values and 

by using the selected states Xi,t to define the region in which pj(x) has extreme values. 

8. (Currently Amended) Method as claimed in one of claims wherein parameters of the 
sampling distribution function pj{x) of at least one iteration are determined by a linear least 
square fit of the logarithm of the un-normalized sampling distribution function pj(x) to the 
logarithm of the product p,tO(xi,t) . 

9. (Currently Amended) Method as claimed in one of claims l-te-8, wherein the normalization 
constant of the sampling distribution function pj(x) of at least one iteration is estimated from 
the sampled states x i { and their weighting factors p,t . 

10. (Currently Amended) Method as claimed in one of claims l-te-9, wherein at least three 
iterations are done. 

1 1. (Currently Amended) Method as claimed in one of claims 1 to 10 , wherein the function 
0(x) is a function of a set of at least two functions Q={Oi(x) s 02(x) ) . . (.*)} , pref e rably 
including at least one of the following functions 

at least one property for which at least one estimate is derived in the analysis of the fifth step, 
at least one function that is large for states that must be sampled to e nsure transitions b e tween 
important regions, 

the inverse of the probability distribution function of at l e ast on e property of the system, and 
the inverse of the probability distribution of the negativ e logarithm of the distribution function 
of the system . 
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12. (Currently Amended) A computer software product for sampling a state space by iteratively 
generating states and their weighting factors fk$ , wherein the index i is the iteration 
parameter and the index / distinguishes different states Xi,t generated by an iteration i, said 
product is characterized bv including a computer-readable medium in which program 
instructions are stored, which instructions, when read by a computer, enable the computer toi 
select an initial sampling distribution function pi(x), 

to execute a fifth step for performing an analysis^ 
and to execute an iteration procedure including, 

a second step for generating Nj states xjj by a numerical sampling algorithm and 
a fourth step for testing at least one criterion to decide whether to continue the iteration 
procedure or to stop the iteration procedure and to go to a fifth step in order to perform the 
analysis using the simulated data, wherein 
characterized in that the iteration further includes 

a third step determining weighting factors pi,t for states Xi,t generated so far by using 
sampling distribution functions p t (x) determined so far and 

a fitting step for determining a sampling distribution function pj(x) for the next iteration by 
fitting pj(x) to pi,tO(xi,t) for states Xi,t generated so far, wherein 0(xis) is a function, 
respectively a property, of the states xi,t . 

13. (New) Method as claimed in claim 2, wherein the sampling distribution function pj(x) of at 
least one iteration is fitted such that the sampling distribution function pj{x) is large for at 
least one state Xi,t with a large product pi,tO{xi.t) , and tends to be small for states with a 
small product p,tO(xi,t) . 

14. (New) Method as claimed in claim 1, wherein the sampling distribution function pj(x) of at 
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least one iteration is a function with at least one constraint, having an extreme value in the 

region of at least one selected state xj.constr and having vanishing values at a distance from the 

at least one selected state Xj&mtr . 

15. (New) Method as claimed in claim 1, wherein the sampling distribution function pj{x) of at 
least one iteration is a function with at least one constraint, having an extreme value in the 
region of at least one selected state xjxonstr and having vanishing values at a distance from the 
at least one selected state xjjnmtr , wherein the constraint is a harmonic constraint with at least 
one constant kconstr defining the strength of the constraint at the selected state xj,constr . 

16. (New) Method as claimed in claim 1, wherein the numerical sampling algorithm of at least 
one iteration generates correlated states xjj , wherein the sampling distribution function 
pj(x) has a maximum that biases sampling into a sub-region of the state space, wherein the 
sampling distribution function pj(x) is fitted such that the maximum is in a region where the 
product pi,tO(xij) has a maximum, and wherein the sampling starts from a state close to the 
maximum of the sampling distribution function pj(x) . 

17. (New) Method as claimed in claim 1, wherein the function 0(x) is a function of a set of at 
least two functions ©={Ch(x) i 02(x) ) . . .,0aw(x)}, including at least one of the following 
functions: 

at least one property for which at least one estimate is derived in the analysis of the fifth step, 
at least one function that is large for states that must be sampled to ensure transitions between 
important regions, 

the inverse of the probability distribution function of at least one property of the system, and 
the inverse of the probability distribution of the negative logarithm of the distribution function 
of the system. 
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18. (New) A computer-readable medium comprising executable program instructions configured 
to cause a computer to perform the method of claim 1. 



19. (New) A computer program, adapted to cause a computer to perform the method of claim 1. 



20. (New) A computer-readable medium comprising the computer program of claim 19. 
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